function [F] = my_DFT(fp)
    % 获取图像的大小
    [M, N] = size(fp);
    
    % 创建频域变量
    F = zeros(M, N);
    
    % 遍历图像的每一个像素点，计算傅里叶变换
    for u = 1:M
        for v = 1:N
            % 计算频域中的每个点
            sum = 0;
            for x = 1:M
                for y = 1:N
                    sum = sum + fp(x, y) * exp(-2*pi*1i * ((u-1)*(x-1)/M + (v-1)*(y-1)/N));
                end
            end
            F(u, v) = sum;
        end
    end
end
